System and method for calibrating light intensity

ABSTRACT

Systems and methods are provided for calibrating light intensity. An exemplary method for light intensity calibration may comprise: obtaining a plurality of intensity distributions of reflected light from an area, wherein each of the intensity distributions is associated with a beam; determining a reference intensity distribution from the plurality of intensity distributions, wherein the reference intensity distribution is associated with a reference beam, the plurality of intensity distributions excluding the reference intensity distribution are non-reference intensity distributions, and the non-reference intensity distributions are each associated with a non-reference beam; and aligning each of the non-reference intensity distributions to the reference distribution to calibrate the non-reference intensity distributions against the reference distribution.

CROSS REFERENCE TO RELATED APPLICATION

This application is a continuation application of U.S. Non-provisionalpatent application Ser. No. 15/886,206, filed Feb. 1, 2018 and titled“SYSTEM AND METHOD FOR CALIBRATING LIGHT INTENSITY,” the entire contentsof which are incorporated herein by reference.

TECHNICAL FIELD

This disclosure generally relates to methods and devices for lightintensity calibration, and in particular, to methods and devices forLIDAR (Light Detection And Ranging) intensity calibration.

BACKGROUND

Sensor devices have been widely used for detection and ranging. Onecritical step for ensuring accurate performances is calibration. Throughcalibration, erroneous intensity readings can be corrected to a standardlevel.

SUMMARY

Various embodiments of the present disclosure can include systems,methods, and non-transitory computer readable media for light intensitycalibration. According to one aspect, an exemplary method for lightintensity calibration may comprise obtaining a plurality of intensitydistributions of reflected light from an area, wherein each of theintensity distributions is associated with a beam; determining areference intensity distribution from the plurality of intensitydistributions, wherein the reference intensity distribution isassociated with a reference beam, the plurality of intensitydistributions excluding the reference intensity distribution arenon-reference intensity distributions, and the non-reference intensitydistributions are each associated with a non-reference beam; andaligning each of the non-reference intensity distributions to thereference distribution to calibrate the non-reference intensitydistributions against the reference distribution. In some embodiments,the method may further comprise obtaining an intensity mapping that mapsany raw non-reference beam intensity reading to a calibrated intensityreading.

In some embodiments, aligning each of the non-reference intensitydistributions to the reference distribution may comprise applying aDynamic Time Warping (DTW) algorithm to align each of the non-referenceintensity distributions to the reference distribution.

In some embodiments, a LIDAR device may comprise N rotatable lightemitters and detectors mountable on a movable vehicle, the each lightemitter emitting the each beam and the each light detector detecting theeach corresponding light intensity while in rotation as the vehiclemoves relative to the area. The beams may be Beams 1, 2, . . . i−1, i,i+1, i+2, . . . , and i+(N−i) in an order based on a physical proximityamong intersections of the beams with the area at the same time, i beingany number from 1 to N. The reference beam may be Beam_(i). Aligningeach of the non-reference intensity distributions to the referencedistribution may comprise, for the intensity distributions, recursivelyaligning the each non-reference beam to a physically closest beamtowards the reference beam until aligning to the reference beam.

In some embodiments, for the intensity distributions, recursivelyaligning the each non-reference beam to the physically closest beamtowards the reference beam may comprise aligning Beam_(i+1) to Beam_(i)and Beam_(i−1) to Beam_(i); aligning Beam_(i+2) to Beam_(i+1) andBeam_(i−2) to Beam_(i−1); repeating the alignments until each of thenon-reference beams is aligned to another beam towards the referencebeam; and propagating the alignments throughout the non-reference beamsto align each of the non-reference beams to the reference beam.

In some embodiments, the physical proximity among the intersections ofthe beams with the area may correspond to another physical proximityamong the light emitters on the LIDAR device.

In some embodiments, the light beams may be respectively emitted from aplurality of light emitters of a LIDAR device. The LIDAR device mayfurther comprise a plurality of light detectors for respectivelyreceiving the reflected light. In association with the each beam, theeach light emitter and associated light detector may be movable relativeto the area and configured to obtain an intensity reading at varioustime points to obtain the corresponding intensity distribution over thearea.

In some embodiments, the LIDAR device may be disposed on a vehicle. Theeach light emitter and associated light detector may be configured torotate relative to the vehicle as the vehicle moves relative to thearea. In association with the each beam, the each intensity distributionover the area may comprise a plurality of the intensity readings causedby the rotation and the vehicle movement.

In some embodiments, each of the intensity distributions may berepresentable in a histogram associating an light intensity reading anda frequency of occurrence of the light intensity reading in the area.The area may comprise a marker corresponding to one or more intensitypeaks in one or more histograms of the intensity distributions. Thereference intensity distribution may be representative of the marker ina background.

According to another aspect, a system for light intensity calibrationmay comprise a processor configured to obtain a plurality of intensitydistributions of reflected light from an area, wherein each of theintensity distributions is associated with a beam; determine a referenceintensity distribution from the plurality of intensity distributions,wherein the reference intensity distribution is associated with areference beam, the plurality of intensity distributions excluding thereference intensity distribution are non-reference intensitydistributions, and the non-reference intensity distributions are eachassociated with a non-reference beam; and align each of thenon-reference intensity distributions to the reference distribution tocalibrate the non-reference intensity distributions against thereference distribution.

According to yet another aspect, the present disclosure provides anon-transitory computer-readable storage medium coupled to a processorand comprising instructions that, when executed by the processor, causethe processor to perform a method for light intensity calibration. Themethod may comprise obtaining a plurality of intensity distributions ofreflected light from an area, wherein each of the intensitydistributions is associated with a beam; determining a referenceintensity distribution from the plurality of intensity distributions,wherein the reference intensity distribution is associated with areference beam, the plurality of intensity distributions excluding thereference intensity distribution are non-reference intensitydistributions, and the non-reference intensity distributions are eachassociated with a non-reference beam; and aligning each of thenon-reference intensity distributions to the reference distribution tocalibrate the non-reference intensity distributions against thereference distribution.

According to yet another aspect, a system for light intensitycalibration, implementable on a LIDAR (Light Detection And Ranging)device, the system may comprise a plurality of rotatable light emittersand detectors mountable on a movable vehicle and a processor coupled tothe detectors. The each light emitter may emit a beam and the each lightdetector may detect an reflection of the corresponding beam while inrotation as the vehicle moves relative to an area. The processor may beconfigured to: obtain a plurality of intensity distributions of thereflected light from the area; determine a reference intensitydistribution from the plurality of intensity distributions, wherein thereference intensity distribution is associated with a reference beam,the plurality of intensity distributions excluding the referenceintensity distribution are non-reference intensity distributions, andthe non-reference intensity distributions are each associated with anon-reference beam; and align each of the non-reference intensitydistributions to the reference distribution to calibrate thenon-reference intensity distributions against the referencedistribution.

In some embodiments, to align the non-reference intensity distributionto the reference distribution, the processor is configured torecursively align the non-reference intensity distribution to aneighboring non-reference intensity distribution towards the referenceintensity distribution until aligned with the reference intensitydistribution. The plurality of light emitters and detectors may compriseN pairs of light emitters and detectors associated with N emitted beamsand N reflections of the emitted beams. The N beams may be Beams 1, 2, .. . i−1, i, i+1, i+2, . . . , and i+(N−i) in an order based on aphysical proximity, i being any number from 1 to N, and the referencebeam being Beam_(i). In one example, the physical proximity may be basedon intersections of the beams with the area at the same time. In anotherexample, the physical proximity may be based on emitters or detectors ofthe beams on the system.

These and other features of the systems, methods, and non-transitorycomputer readable media disclosed herein, as well as the methods ofoperation and functions of the related elements of structure and thecombination of parts and economies of manufacture, will become moreapparent upon consideration of the following description and theappended claims with reference to the accompanying drawings, all ofwhich form a part of this specification, wherein like reference numeralsdesignate corresponding parts in the various figures. It is to beexpressly understood, however, that the drawings are for purposes ofillustration and description only and are not intended as a definitionof the limits of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

Certain features of various embodiments of the present technology areset forth with particularity in the appended claims. A betterunderstanding of the features and advantages of the technology will beobtained by reference to the following detailed description that setsforth illustrative embodiments, in which the principles of the inventionare utilized, and the accompanying drawings of which:

FIG. 1 illustrates an exemplary system for light intensity calibration,in accordance with various embodiments.

FIG. 2 illustrates another exemplary system for light intensitycalibration, in accordance with various embodiments.

FIG. 3A illustrates an exemplary system for LIDAR intensity calibration,in accordance with various embodiments.

FIG. 3B illustrates an exemplary LIDAR system, in accordance withvarious embodiments.

FIG. 3C illustrates exemplary intensity spots for multiple rotating andstationary beams, in accordance with various embodiments.

FIG. 3D illustrates exemplary intensity spots for a rotating and movingbeam, in accordance with various embodiments.

FIG. 3E illustrates exemplary intensity distributions and calibration,in accordance with various embodiments.

FIG. 3F illustrates an exemplary intensity mapping, in accordance withvarious embodiments.

FIG. 3G illustrates three exemplary sets of raw intensity output andcalibrated intensity output comparisons, in accordance with variousembodiments.

FIG. 4 illustrates a flowchart of an exemplary method for lightintensity calibration, in accordance with various embodiments.

FIG. 5 illustrates a flowchart of an exemplary method for lightintensity alignment, in accordance with various embodiments.

FIG. 6 illustrates a block diagram of an exemplary computer system inwhich any of the embodiments described herein may be implemented.

DETAILED DESCRIPTION

Multiple-beam LIDARs are critical sensors responsible for localization,mapping, and various other functions of autonomous vehicles. Anexemplary LIDAR system may comprise multiple sets of light emitters(e.g., lasers) and light detectors, each set configured to emit a beamtowards a distant object and receive reflected light signals. Based onthe received signals, the LIDAR can obtain (x, y, z) coordinates andintensity values corresponding to the distant object.

In some implementations, the distant object may comprise traffic signsor other important information, and accurate LIDAR detection isessential for correctly recognizing such information. To that end, theintensity values need to be calibrated across the beams. For example,calibrated beams should correctly give the same intensity reading on thesame target (e.g., a part of road surface). In current technologies, oneapproach to the problem divides the target space into small “cells,” andfor any beam r that returns value i on cell p, uses the consensus ofother beams in p as an estimation. This existing method reliessignificantly on accurate pose estimations of the LIDAR in order toregister detection points from multiple beams correctly into the samecell.

The disclosed method obviates the requirement for accurate poseestimation by aligning intensity distribution over a much larger area(e.g., on the order of hundreds of meters) where all the beams share asignificant overlap, while achieving an accurate calibration. Thedisclosed methods are also unsupervised and convenient forimplementation.

FIG. 1 illustrates an exemplary system 100 for light intensitycalibration, in accordance with various embodiments. As shown in FIG. 1,the exemplary system 100 can comprise at least one computing system 102that includes one or more processors 104 and memory 106. The memory 106may be non-transitory and computer-readable. The memory 106 may storeinstructions that, when executed by the one or more processors 104,cause the one or more processors 104 to perform various operationsdescribed herein. The system 102 may be implemented on or as variousdevices such as mobile phone, tablet, server, computer, wearable device(smart watch), etc. The system 102 may be installed with appropriatesoftware (e.g., data transfer program, etc.) and/or hardware (e.g., wireconnections, wireless connections, etc.) to access other devices of thesystem 100.

The system 100 may include one or more data stores (e.g., a data store108) and one or more computing devices (e.g., a computing device 109)that are accessible to the system 102. In some embodiments, the system102 may be configured to obtain data (e.g., location, time, movingspeed, acceleration, map) from the data store 108 (e.g., a clouddatabase) and/or the computing device 109 (e.g., a server, a mobilephone, a vehicle computer).

The system 100 may further include a system 101 (e.g., a LIDAR system)coupled to the system 102. The system 101 may comprise one or more lightdetectors (e.g., light detectors 110 a, 111 a, etc.) and one or morelight emitters (e.g., light emitters 110 b, 111 b, etc.). The lightemitters and detectors (e.g., 64-channel infrared laser and detector)may form emitter-detector pairs. Alternatively, one emitter maycorrespond to multiple detectors and vice versa. In some embodiments,the system 102 (e.g., the processor 104) may be configured to obtaindata (e.g., receive intensity readings, 3D locations of the reflectionspots, time points of the intensity readings) from the one or more lightdetectors, and perform various light intensity calibration methods andsteps described herein. Each intensity reading may correspond to a beamand a light emitter that emitted the beam. The calibration method mayyield an intensity mapping for each detector-emitter pair detection beamto map any raw intensity reading from this beam to a calibratedintensity reading.

In some embodiments, the system 102 and the one or more of the computingdevices (e.g., the computing device 109) may be integrated in a singledevice or system. Alternatively, the system 102 and the one or morecomputing devices may operate as separate devices. The one or more datastore (e.g., the data store 108) may be anywhere accessible to thesystem 102, for example, in the memory 106, in the computing device 109,in another device (e.g., network storage device) coupled to the system102, or another storage location (e.g., cloud-based storage system,network file system, etc.), etc. The system 102 may be implemented as asingle system or multiple systems coupled to each other. In general, thesystem 102, the computing device 109, the data store 108, and the system101 may be able to communicate with one another through one or morewired or wireless networks (e.g., the Internet) through which data canbe communicated.

FIG. 2 illustrates an exemplary system 200 for light intensitycalibration, in accordance with various embodiments. The operationsshown in FIG. 2 and presented below are intended to be illustrative. Thesystem 200 is similar to the system 100 described above, except that thesystem 102 and the system 101 illustrated in FIG. 1 are merged to obtainthe system 103 in FIG. 2. The system 103 may be implemented as a LIDARsystem. The system 103 may comprise a processor 104, a memory 106, oneor more light detectors (e.g., 110 a, 111 a, etc.), and one or morelight emitters (e.g., 110 b, 111 b, etc.), each configured to performsteps similar to those described above with reference to FIG. 1. Thedata store 108 and the computing device 109 may be coupled or integratedto the system 103. Various aspects of the system 100 and the system 200are described below in reference to FIG. 3A to FIG. 6.

FIG. 3A illustrates an exemplary system for LIDAR intensity calibration,in accordance with various embodiments. The operations shown in FIG. 3Aand presented below are intended to be illustrative. As shown in theside view of FIG. 3A, a system 312 may be disposed on a vehicle 311(e.g., car, taxi, internet-based service vehicle, motorcycle, bus, boat,drone, etc.). The system 312 may comprise the system 101 or the system103 described above. If the system 312 comprises the system 101, thesystem 102 may be disposed also in the system 312, disposed in thevehicle 311, or remotely (e.g., wirelessly) coupled to the system 312.The vehicle 311 may be, but not necessarily, autonomous. That is, thevehicle 311 may be capable of navigating without a human driver. Thevehicle 311 may be configured to move around similar to regular cars.The driverless navigation may be achieved at least based in part on thesensing output by the system 312. The system 312 may be configured torotate (e.g., with respect to an axis vertical to the vehicle 311).

In various embodiments, the system 312 may emit one or more beams (e.g.,Beam_(N) to Beam₁ from various light emitters) towards a distant targetand correspondingly detect the reflected beams. The emitted beams maytouch upon and reflect off any target, such as a distant area 318 on theground, an area 319 on an erected board, etc. For example, Beam₁ mayreflect off the area 318 at Intersection₁, Beam₂ may reflect off thearea 318 at Intersection₂, and Beam N may reflect off the area 318 atIntersection_(N). The physical proximity between the intersections maycorrespond to the physical proximity between the beams. For example,Intersection₁ is closer to Intersection₂ than Intersection_(N), whichcorresponds to the closer proximity between Beam₁ and Beam₂. Further,the area scanned by the beams may comprise one or more signs, markings,labels, or the equivalents that carry information (e.g., traffic signs,map information). By the scan, intensity readings off the spots on thesigns and the surrounding background can be obtained, and intensitydifferences can be used to extract the carried information.

FIG. 3B illustrates an exemplary LIDAR system, in accordance withvarious embodiments. The operations shown in FIG. 3B and presented beloware intended to be illustrative. FIG. 3B illustrates further details ofthe system 312 according to some embodiments in side view. The system312 may comprise multiple light emitters and light detectors configuredto detect reflections from multiple beams at the same time. In someembodiments, the light emitters and detectors may be paired (e.g., lightemitter-detector₁ to light emitter-detector_(N)), each emitting anoutgoing beam and receiving an incoming beam as a reflection of theoutgoing beam. Although in this figure a light emitter and lightdetector pair is shown as an integrated component, an independent lightemitter and an independent light detector can be disposed as separatecomponents and pairing for detection. Further, one emitter may bedetected by multiple detectors, and one detector may detect multiplereflections from multiple beams. As shown in this figure, the physicalproximity between the beams may correspond to the physical proximitybetween the corresponding light emitter-detectors on the system 312. Forexample, the light emitter-detector pairs may be arranged on a sidesurface from top to bottom (from N to 1) facing towards a distantground. Thus, in view of FIG. 3A and FIG. 3B, the physical proximityamong the intersections of the beams with the area 318 may correspond toanother physical proximity among the light emitters on the LIDAR device.

FIG. 3C illustrates exemplary intensity spots for multiple rotating andstationary beams, in accordance with various embodiments. The operationsshown in FIG. 3C and presented below are intended to be illustrative. Inthe top view of the FIG. 3C, the vehicle 311 is relatively stationary toa distant area 338 on the ground, which comprises a “stop” sign. In someembodiments, the system 312 may be configured to rotate with respect tothe vehicle 311, while the light emitter-detector pairs respectivelyemit light and capture reflected light intensities. By the rotation,each light beam facing the ground scans a circumference of a circlecentered at the system 312 (e.g., Ring₁ to Ring_(N)). Since the lightemitters may be pulsed lasers, intensities of reflections from variousspots along the circumference are measured. By increasing the spotdensity, reflections off nearly the entire circumference can bemeasured. Further, the light emitter-detectors can be positioned atvarious angles with respect to the system 312 such that reflection spotsat various distances from the vehicle 311 can be measured. For example,in this figure, Ring₁ comprises reflection spots closer to the vehicle311 than Ring_(N). The distances among the rings can be adjusted.

FIG. 3D illustrates exemplary intensity spots for a rotating and movingbeam, in accordance with various embodiments. The operations shown inFIG. 3D and presented below are intended to be illustrative. In the topview of the FIG. 3D, the vehicle 311 may move relatively to (e.g.,towards) the distant area 338. Although multiple beams may be emittedsimultaneously from the system 312 to detect their reflectedintensities, only one beam is shown in this figure to keep theillustration clean. Here, the system 312 may be configured to performthe scan described above with reference to FIG. 3C continuously as thevehicle 311 moves towards the area 338 (e.g., from vehicle 311 a andsystem 312 a at time t₁ to vehicle 311 b and system 312 b at time t₂).At time t₁, an exemplary Beam₁ may rotate and yield Ring₁ of t₁, a partof which barely enters a closer edge of the area 338. At time t₂, thesame Beam₁ may rotate and yield Ring₁ of t₂, a part of which is leavinga further edge of the area 338. Between the time t₁ and t₂ (not shownhere), the Beam₁ may have been continuously scanning, and some of thedetected intensity are reflected off spots across the entire area 338.Similar scanning and intensity detections can be obtained from otherbeams of the system 312.

In some embodiments, as shown above, the scanning across the area 338can be enabled by the rotation and movement of the vehicle. The lightbeams may be fixed to the vehicle 311 except for the rotation.Nevertheless, there can be many other mechanisms for alternativelyachieving the scanning across the area. In various examples, the lightemitters may be configured to move translationally relative to thevehicle, the light emitters may be configured to rotate in a planevertical to the vehicle, the emitters may comprise optics configured totune the angle of the emitted beam relative to vehicle, etc.

FIG. 3E illustrates exemplary intensity distributions and calibration,in accordance with various embodiments. The operations shown in FIG. 3Eand presented below are intended to be illustrative. Continuing fromFIG. 3D, reflection spots of Beam₁ in the area 338 may be aggregatedbetween t₁ and t₂ as the system 312 rotates and the vehicle 311advances. The aggregated spot distribution for Beam₁ in the area 338 isshown in Diagram (1) of FIG. 3E. As discussed above, since eachintensity measurement may be associated with a time of the measurementand/or a 3D-location of the reflection spot, reflection spots that fallin a certain area (e.g., the area 338 or some designated portion of thearea 338) may be determined based on the 3D-locations and a boundary ofthe area, or based on the time stamps (e.g., t₁ and t₂). In practice,the reflection spots can be much denser over an area to increaseaccuracy.

Once the reflection spots within the area 338 are determined, as shownin Diagram (2) of FIG. 3E, an intensity distribution can be obtained forBeam₁. The intensity distribution can be represented, for example, by ahistogram with the intensity on the x-axis and the frequency on they-axis (here the histogram is smoothed to a curve). That is, each of theintensity distributions may be representable in a histogram associatingan light intensity reading and a frequency of occurrence of the lightintensity reading in the area. The intensity axis may, for example,correspond to a pixel detection range from 0 (black) to 255 (white). Thearea 338 may comprise a marker (e.g., the “stop”) corresponding to oneor more intensity peaks in one or more histograms of the intensitydistributions. In one exemplary embodiment, the reference intensitydistribution may be representative of the area, with some of theintensity peaks representing reflections off the marker and some otherintensity peaks representing reflections off the surrounding background.For example, the reflection intensities off the “stop” sign maycorrespond to the left peak, and reflection intensities off thesurrounding background of the “stop” sign may correspond to the rightpeak.

Similarly, as shown in Diagram (3) of FIG. 3E, intensity distributionsfor all beams from Beam₁ to Beam_(N) can be obtained. That is, aplurality of intensity distributions of reflected light from the area338 can be obtained, wherein each of the intensity distributions isassociated with a beam. As shown in this diagram, points along thehistograms of various beams may be deformed relative to each other, eventhough the reflection spots for each beam covering the same area 338 aresupposed to yield similar intensity histograms. The deformation may ormay not be uniform. That is, the deformation may appear random acrossthe histograms. A relative deformation between two neighboringhistograms may comprise shifting in the x and/or y directions,stretching or squeezing, mismatch in peak heights, or any otherdeformation obtained by the device measurement. In various embodiments,the histograms of beams may be more complicated than those shown in thefigure, and the histograms may have larger differences. For example, theshifts from Beam₁ to Beam_(N) may appear random. Beam₂ may appearshifted left from Beam₁, Beam₃ may appear shifted right from Beam₂, andBeam₃ may have one peak (e.g., its right peak) missing (e.g., below athreshold such that it is noticeable). The corresponding peaks in thehistograms may have different scalers in one direction (e.g., fromBeam_(i) to Beam₁, from Beam_(i) to Beam_(N)). As such, the calibrationmay be non-linear. Thus, each of these intensity distributions need tobe calibrated such that the intensity measurements of different beamsacross the same area yield similar (e.g., completely overlapping)intensity distributions. For clearer illustration, the following Diagram(4) of FIG. 3E may only show shifted histograms as an example.

With reference to Diagram (4) of FIG. 3E, a reference intensitydistribution can be determined from the plurality of intensitydistributions, wherein the reference intensity distribution isassociated with a reference beam, the plurality of intensitydistributions excluding the reference intensity distribution arenon-reference intensity distributions, and the non-reference intensitydistributions are each associated with a non-reference beam. Each of thenon-reference intensity distributions can be aligned to the referencedistribution to calibrate the non-reference intensity distributionsagainst the reference distribution. There can be many ways to determinethe reference intensity distribution. For example, if the area consistsof a uniform sign with a uniform background, the reference intensitydistribution should comprise two non-overlapping peaks representingreflected intensities off the two distinct regions, and an intensitydistribution having such well-positioned characteristics can bedetermined as the reference intensity distribution. Exemplary rules suchas peak separation thresholding, peak value range thresholding, and peakheight thresholding can be applied to make such determination. Foranother example, if a predetermined intensity distribution or statisticsof a predetermined intensity distribution for the same area isaccessible to the system 312 (e.g., stored in the system 312 or remotelyavailable), the intensity distribution may be determined as thereference intensity distribution. Exemplary rules such as theleast-mean-square difference can be applied to make such determination.

In some embodiments, as described above, the system 312 (e.g., a LIDARdevice) may comprise N rotatable light emitters and detectors mountableon a movable vehicle, the each light emitter emitting the each beam andthe each light detector detecting the each corresponding lightintensity, while the emitters and detectors are in rotation and thevehicle moves relative to the area. The beams may be Beams 1, 2, . . .i−1, i, i+1, i+2, . . . , and i+(N−i) in an order based on a physicalproximity (e.g., a physical proximity among intersections of the beamswith the area at the same time, a physical proximity among the emittedbeams in the air, a physical proximity among the emitted beams in theair, a physical proximity among the reflected beams in the air, aphysical proximity among the emitters on the LIDAR, a physical proximityamong the detectors on the LIDAR, etc.), i being any number from 1 to N.The reference beam may be beam i.

For example, referring to FIG. 3C, Ring₁ is closer to Ring₂ thanRing_(N). Referring back to FIG. 3E, in some embodiments, aligning eachof the non-reference intensity distributions to the referencedistribution may comprise, for the intensity distributions, recursivelyaligning the each non-reference beam to a physically closest beamtowards the reference beam until aligning to the reference beam. Here,the alignment of beams may refer to the alignment of the intensitydistributions corresponding to the beams. For example, as shown inDiagram (4) of FIG. 3E, Beam_(i+1) may be aligned to Beam_(i), andBeam_(i−1) may be aligned to Beam_(i). For example, a Dynamic TimeWarping (DTW) algorithm can be applied to align the intensitydistributions. Some of the mapping relations are shown as dash arrows inDiagram (4). The DTW algorithm can calculate an optimal match betweentwo given sequences (e.g., in a time series). The sequences are “warped”non-linearly in one dimension to determine a measure of their similarityindependent of certain non-linear variations in the dimension.Notwithstanding, a person of ordinary skill in the art should appreciatethat various other alternative methods or algorithms can be applied toalign the intensity distributions.

Further, Beam_(i+2) may be aligned to Beam_(i+1), and Beam_(i−2) may bealigned to Beam_(i−1). The alignments may be repeated until each of thenon-reference beams is aligned to another beam towards the referencebeam (e.g., a next beam towards the reference beam). Further, thealignments can be propagated throughout the non-reference beams to aligneach of the non-reference beams to the reference beam. For example,given that Beam_(i+1) has been aligned to Beam_(i) and that Beam_(i+2)has been aligned to Beam_(i+1), the propagation can pass the alignmentof Beam_(i+1) to Beam_(i+2), such that Beam_(i+2) is aligned toBeam_(i). Similarly, any alignment of a non-reference beam can be passedon and used in aligning one or more other non-reference beams furtheraway from the reference beam. Aligning neighboring beams and propagatingthe alignment can help increase the alignment accuracy and decrease thecomputation power consumption, because neighboring beams cover closerareas (due to the physical proximity described above) and deviationsbetween the neighboring beams are smaller than those between thenon-neighboring beams. Thus, by one alignment for the each non-referencebeam and the propagation, all non-reference beams can be aligned to thereference beam. And by aligning all non-reference intensitydistributions, any value on a non-reference histogram can be mapped to avalue (e.g., a calibrated value) on the reference histogram.

In some embodiments, corresponding transformation steps can begeneralized according to the alignment described above. For example, forBeam_(i−1) (a non-reference beam) and Beam_(i) (a reference beam),alignment function W(k; i−1 to i) can be obtained for Beam_(i−1), whichconvert Beam_(i−1)'s raw intensity k to intensity W(k; i−1 to i) inalignment with Beam i. In one example, Beam₁₅ is a reference beam, andits intensity distribution is a reference to which the intensitydistributions of other beams are to be aligned. If a frequency value ofa peak in Beam₁₅ is 100, and Beam₁₄ has a similar peak around a similarintensity value with a frequency value k of 90, alignment function W(90;14 to 15) would yield a calibrated intensity of 100 for the peak inBeam₁₄. A similar alignment function W(k; i+1 to i) can be used forbeams higher in number than the reference beam. For example, if afrequency value of a peak in Beam₁₅ is 100, and Beam₁₆ has a similarpeak around a similar intensity value with a frequency value k of 110,the alignment function W(110; 16 to 15) would yield an calibratedintensity of 100 for the peak in Beam₁₆.

Further, in the calibration propagation step, the global calibrationfunction can be obtained by iterating the alignment function over allintermediate beams. Continuing the example of alignment function W(90;14 to 15) described above, to align a similar peak in Beam₁₃, the globalcalibration function this alignment is W(k; 13 to 15). This calibrationcomprises iterations of the alignment over Beam₁₄ and Beam₁₅, that is,W(k; 13 to 15)=W(k; 13 to 14) in combination with W(k; 14 to 15). Thus,W(k; 13 to 15) can be obtained from W(W(k; 13 to 14); 14 to 15).Similarly, any of the non-reference beam's intensity value can becalibrated by iterating the k value over intermediate beams between thetarget bream and the reference beam.

FIG. 3F illustrates an exemplary intensity mapping, in accordance withvarious embodiments. The operations shown in FIG. 3F and presented beloware intended to be illustrative. Continuing from FIG. 3E, an intensitymapping that maps any raw non-reference beam intensity reading to acalibrated intensity reading can be obtained based on the intensitydistribution alignment for each non-reference beam. For eachnon-reference beam, various intensity values may correspond to the sameor different alignment functions. In FIG. 3F, three exemplary alignmentfunctions are plotted for three raw intensity ranges. Once this mappingis obtained, any future raw intensity reading on the x-axis can bemapped to a calibrated intensity reading on the y-axis. Although in thisfigure the intensity mapping is shown as a 2D plot, various otherintensity mapping representations can be alternatively used.

FIG. 3G illustrates three exemplary sets of raw intensity output andcalibrated intensity output comparisons, in accordance with variousembodiments. The operations shown in FIG. 3G and presented below areintended to be illustrative. The left image of each set shows the rawintensity output without applying the disclosed methods, and the rightimage of each set shows the corresponding intensity output afterapplying the disclosed methods. The intensity outputs in the rightimages all have shown greater clearness, smoothness, and accuracy atleast with respect to road imagery (e.g., lane marker, direction, sign)represented as white lines. Thus, the disclosed systems and methods canbe used to calibrate LIDARs for better performance on roads.

FIG. 4 illustrates a flowchart of an exemplary method 400, according tovarious embodiments of the present disclosure. The method 400 may beimplemented in various environments including, for example, the system100 of FIG. 1 or the system 200 of FIG. 2. The exemplary method 400 maybe implemented by one or more components of the system 102 or 103 (e.g.,the processor 104, the memory 106). The exemplary method 400 may beimplemented by multiple systems similar to the system 102 or 103. Theoperations of method 400 presented below are intended to beillustrative. Depending on the implementation, the exemplary method 400may include additional, fewer, or alternative steps performed in variousorders or in parallel.

At block 402, a plurality of intensity distributions of reflected lightfrom an area may be obtained, wherein each of the intensitydistributions is associated with a beam. At block 404, a referenceintensity distribution from the plurality of intensity distributions maybe determined, wherein the reference intensity distribution isassociated with a reference beam, the plurality of intensitydistributions excluding the reference intensity distribution arenon-reference intensity distributions, and the non-reference intensitydistributions are each associated with a non-reference beam. At block406, each of the non-reference intensity distributions may be aligned tothe reference distribution to calibrate the non-reference intensitydistributions against the reference distribution. To align anynon-reference intensity distribution to the reference distribution, thenon-reference intensity distribution may be recursively aligned to aneighboring non-reference intensity distribution towards the referenceintensity distribution until aligned with the reference intensitydistribution. In some embodiments, a LIDAR (Light Detection And Ranging)device may comprise N rotatable light emitters and N detectors mountableon a movable vehicle, the each light emitter emitting the each beam andthe each light detector detecting the each corresponding light intensitywhile in rotation as the vehicle moves relative to the area. The beamsmay be Beams 1, 2, . . . i−1, i, i+1, i+2, . . . , and i+(N−i) in anorder based on a physical proximity among intersections of the beamswith the area at the same time, i being any number from 1 to N. Thereference beam may be Beam_(i). Aligning each of the non-referenceintensity distributions to the reference distribution may comprise, forthe intensity distributions, recursively aligning the each non-referencebeam to a physically closest beam towards the reference beam untilaligning to the reference beam. For the intensity distributions, asshown in FIG. 5, recursively aligning the each non-reference beam to thephysically closest beam towards the reference beam may comprise: (1)aligning Beam_(i+1) to Beam_(i) and Beam_(i−1) to Beam_(i), as shown instep 422, (2) aligning Beam_(i+2) to Beam_(i+1) and Beam_(i−2) toBeam_(i−1), as shown in step 424, (3) repeating the alignments untileach of the non-reference beams is aligned to another beam towards thereference beam, as shown in step 426, and (4) propagating the alignmentsthroughout the non-reference beams to align each of the non-referencebeams to the reference beam, as shown in step 428. Optionally, anintensity mapping that maps any raw non-reference beam intensity readingto a calibrated intensity reading may be further obtained.

In some embodiments, the light beams may be respectively emitted from aplurality of light emitters of a LIDAR device. The LIDAR device mayfurther comprise a plurality of light detectors for respectivelyreceiving the reflected light. In association with the each beam, theeach light emitter and associated light detector may be movable relativeto the area and configured to obtain an intensity reading at varioustime points to obtain the corresponding intensity distribution over thearea.

In some embodiments, the LIDAR device may be disposed on a vehicle. Theeach light emitter and associated light detector may be configured torotate relative to the vehicle as the vehicle moves relative to thearea. In association with the each beam, the each intensity distributionover the area may comprise a plurality of the intensity readings causedby the rotation and the vehicle movement.

The techniques described herein are implemented by one or morespecial-purpose computing devices. The special-purpose computing devicesmay be hard-wired to perform the techniques, or may include circuitry ordigital electronic devices such as one or more application-specificintegrated circuits (ASICs) or field programmable gate arrays (FPGAs)that are persistently programmed to perform the techniques, or mayinclude one or more hardware processors programmed to perform thetechniques pursuant to program instructions in firmware, memory, otherstorage, or a combination. Such special-purpose computing devices mayalso combine custom hard-wired logic, ASICs, or FPGAs with customprogramming to accomplish the techniques. The special-purpose computingdevices may be desktop computer systems, server computer systems,portable computer systems, handheld devices, networking devices or anyother device or combination of devices that incorporate hard-wiredand/or program logic to implement the techniques. Computing device(s)are generally controlled and coordinated by operating system software.Conventional operating systems control and schedule computer processesfor execution, perform memory management, provide file system,networking, I/O services, and provide a user interface functionality,such as a graphical user interface (“GUI”), among other things.

FIG. 6 is a block diagram that illustrates a computer system 500 uponwhich any of the embodiments described herein may be implemented. Thesystem 500 may correspond to the system 102 or 103 described above. Thecomputer system 500 includes a bus 502 or other communication mechanismfor communicating information, one or more hardware processors 504coupled with bus 502 for processing information. Hardware processor(s)504 may be, for example, one or more general purpose microprocessors.The processor(s) 504 may correspond to the processor 104 describedabove.

The computer system 500 also includes a main memory 506, such as arandom access memory (RAM), cache and/or other dynamic storage devices,coupled to bus 502 for storing information and instructions to beexecuted by processor 504. Main memory 506 also may be used for storingtemporary variables or other intermediate information during executionof instructions to be executed by processor 504. Such instructions, whenstored in storage media accessible to processor 504, render computersystem 500 into a special-purpose machine that is customized to performthe operations specified in the instructions. The computer system 500further includes a read only memory (ROM) 508 or other static storagedevice coupled to bus 502 for storing static information andinstructions for processor 504. A storage device 510, such as a magneticdisk, optical disk, or USB thumb drive (Flash drive), etc., is providedand coupled to bus 502 for storing information and instructions. Themain memory 506, the ROM 508, and/or the storage 510 may correspond tothe memory 106 described above.

The computer system 500 may implement the techniques described hereinusing customized hard-wired logic, one or more ASICs or FPGAs, firmwareand/or program logic which in combination with the computer systemcauses or programs computer system 500 to be a special-purpose machine.According to one embodiment, the techniques herein are performed bycomputer system 500 in response to processor(s) 504 executing one ormore sequences of one or more instructions contained in main memory 506.Such instructions may be read into main memory 506 from another storagemedium, such as storage device 510. Execution of the sequences ofinstructions contained in main memory 506 causes processor(s) 504 toperform the process steps described herein. In alternative embodiments,hard-wired circuitry may be used in place of or in combination withsoftware instructions.

The main memory 506, the ROM 508, and/or the storage 510 may includenon-transitory storage media. The term “non-transitory media,” andsimilar terms, as used herein refers to any media that store data and/orinstructions that cause a machine to operate in a specific fashion. Suchnon-transitory media may comprise non-volatile media and/or volatilemedia. Non-volatile media includes, for example, optical or magneticdisks, such as storage device 510. Volatile media includes dynamicmemory, such as main memory 506. Common forms of non-transitory mediainclude, for example, a floppy disk, a flexible disk, hard disk, solidstate drive, magnetic tape, or any other magnetic data storage medium, aCD-ROM, any other optical data storage medium, any physical medium withpatterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, NVRAM, anyother memory chip or cartridge, and networked versions of the same.

The computer system 500 also includes a network interface 518 coupled tobus 502. Network interface 518 provides a two-way data communicationcoupling to one or more network links that are connected to one or morelocal networks. For example, network interface 518 may be an integratedservices digital network (ISDN) card, cable modem, satellite modem, or amodem to provide a data communication connection to a corresponding typeof telephone line. As another example, network interface 518 may be alocal area network (LAN) card to provide a data communication connectionto a compatible LAN (or WAN component to communicated with a WAN).Wireless links may also be implemented. In any such implementation,network interface 518 sends and receives electrical, electromagnetic oroptical signals that carry digital data streams representing varioustypes of information.

The computer system 500 can send messages and receive data, includingprogram code, through the network(s), network link and network interface518. In the Internet example, a server might transmit a requested codefor an application program through the Internet, the ISP, the localnetwork and the network interface 518.

The received code may be executed by processor 504 as it is received,and/or stored in storage device 510, or other non-volatile storage forlater execution.

Each of the processes, methods, and algorithms described in thepreceding sections may be embodied in, and fully or partially automatedby, code modules executed by one or more computer systems or computerprocessors comprising computer hardware. The processes and algorithmsmay be implemented partially or wholly in application-specificcircuitry.

The various features and processes described above may be usedindependently of one another, or may be combined in various ways. Allpossible combinations and sub-combinations are intended to fall withinthe scope of this disclosure. In addition, certain method or processblocks may be omitted in some implementations. The methods and processesdescribed herein are also not limited to any particular sequence, andthe blocks or states relating thereto can be performed in othersequences that are appropriate. For example, described blocks or statesmay be performed in an order other than that specifically disclosed, ormultiple blocks or states may be combined in a single block or state.The exemplary blocks or states may be performed in serial, in parallel,or in some other manner. Blocks or states may be added to or removedfrom the disclosed exemplary embodiments. The exemplary systems andcomponents described herein may be configured differently thandescribed. For example, elements may be added to, removed from, orrearranged compared to the disclosed exemplary embodiments.

The various operations of exemplary methods described herein may beperformed, at least partially, by an algorithm. The algorithm may becomprised in program codes or instructions stored in a memory (e.g., anon-transitory computer-readable storage medium described above). Suchalgorithm may comprise a machine learning algorithm. In someembodiments, a machine learning algorithm may not explicitly programcomputers to perform a function, but can learn from training data tomake a predictions model that performs the function.

The various operations of exemplary methods described herein may beperformed, at least partially, by one or more processors that aretemporarily configured (e.g., by software) or permanently configured toperform the relevant operations. Whether temporarily or permanentlyconfigured, such processors may constitute processor-implemented enginesthat operate to perform one or more operations or functions describedherein.

Similarly, the methods described herein may be at least partiallyprocessor-implemented, with a particular processor or processors beingan example of hardware. For example, at least some of the operations ofa method may be performed by one or more processors orprocessor-implemented engines. Moreover, the one or more processors mayalso operate to support performance of the relevant operations in a“cloud computing” environment or as a “software as a service” (SaaS).For example, at least some of the operations may be performed by a groupof computers (as examples of machines including processors), with theseoperations being accessible via a network (e.g., the Internet) and viaone or more appropriate interfaces (e.g., an Application ProgramInterface (API)).

The performance of certain of the operations may be distributed amongthe processors, not only residing within a single machine, but deployedacross a number of machines. In some exemplary embodiments, theprocessors or processor-implemented engines may be located in a singlegeographic location (e.g., within a home environment, an officeenvironment, or a server farm). In other exemplary embodiments, theprocessors or processor-implemented engines may be distributed across anumber of geographic locations.

Throughout this specification, plural instances may implementcomponents, operations, or structures described as a single instance.Although individual operations of one or more methods are illustratedand described as separate operations, one or more of the individualoperations may be performed concurrently, and nothing requires that theoperations be performed in the order illustrated. Structures andfunctionality presented as separate components in exemplaryconfigurations may be implemented as a combined structure or component.Similarly, structures and functionality presented as a single componentmay be implemented as separate components. These and other variations,modifications, additions, and improvements fall within the scope of thesubject matter herein.

Although an overview of the subject matter has been described withreference to specific exemplary embodiments, various modifications andchanges may be made to these embodiments without departing from thebroader scope of embodiments of the present disclosure. Such embodimentsof the subject matter may be referred to herein, individually orcollectively, by the term “invention” merely for convenience and withoutintending to voluntarily limit the scope of this application to anysingle disclosure or concept if more than one is, in fact, disclosed.

The embodiments illustrated herein are described in sufficient detail toenable those skilled in the art to practice the teachings disclosed.Other embodiments may be used and derived therefrom, such thatstructural and logical substitutions and changes may be made withoutdeparting from the scope of this disclosure. The Detailed Description,therefore, is not to be taken in a limiting sense, and the scope ofvarious embodiments is defined only by the appended claims, along withthe full range of equivalents to which such claims are entitled.

Any process descriptions, elements, or blocks in the flow diagramsdescribed herein and/or depicted in the attached figures should beunderstood as potentially representing modules, segments, or portions ofcode which include one or more executable instructions for implementingspecific logical functions or steps in the process. Alternateimplementations are included within the scope of the embodimentsdescribed herein in which elements or functions may be deleted, executedout of order from that shown or discussed, including substantiallyconcurrently or in reverse order, depending on the functionalityinvolved, as would be understood by those skilled in the art.

As used herein, the term “or” may be construed in either an inclusive orexclusive sense. Moreover, plural instances may be provided forresources, operations, or structures described herein as a singleinstance. Additionally, boundaries between various resources,operations, engines, and data stores are somewhat arbitrary, andparticular operations are illustrated in a context of specificillustrative configurations. Other allocations of functionality areenvisioned and may fall within a scope of various embodiments of thepresent disclosure. In general, structures and functionality presentedas separate resources in the exemplary configurations may be implementedas a combined structure or resource. Similarly, structures andfunctionality presented as a single resource may be implemented asseparate resources. These and other variations, modifications,additions, and improvements fall within a scope of embodiments of thepresent disclosure as represented by the appended claims. Thespecification and drawings are, accordingly, to be regarded in anillustrative rather than a restrictive sense.

Conditional language, such as, among others, “can,” “could,” “might,” or“may,” unless specifically stated otherwise, or otherwise understoodwithin the context as used, is generally intended to convey that certainembodiments include, while other embodiments do not include, certainfeatures, elements and/or steps. Thus, such conditional language is notgenerally intended to imply that features, elements and/or steps are inany way required for one or more embodiments or that one or moreembodiments necessarily include logic for deciding, with or without userinput or prompting, whether these features, elements and/or steps areincluded or are to be performed in any particular embodiment.

The invention claimed is:
 1. A method for light intensity calibration,comprising: obtaining a plurality of intensity distributions ofreflected light from an area, wherein each of the intensitydistributions is associated with a beam; determining a referenceintensity distribution from the plurality of intensity distributions,wherein the reference intensity distribution is associated with areference beam, the plurality of intensity distributions excluding thereference intensity distribution are non-reference intensitydistributions, and the non-reference intensity distributions are eachassociated with a non-reference beam; and aligning each of thenon-reference intensity distributions to the reference distribution byapplying a Dynamic Time Warping (DTW) algorithm to calibrate thenon-reference intensity distributions against the referencedistribution, wherein: a LIDAR (Light Detection And Ranging) devicecomprises N rotatable light emitters and N detectors mountable on amovable vehicle, the each light emitter emitting the each beam and theeach light detector detecting the each corresponding light intensitywhile in rotation as the vehicle moves relative to the area; the beamsare Beams 1, 2, . . . i−1, i, i+1, i+2, . . . , and i+(N−i) in an orderbased on a physical proximity among intersections of the beams with thearea at the same time, i being any number from 1 to N; the referencebeam is Beam_(i); aligning each of the non-reference intensitydistributions to the reference distribution comprises, for the intensitydistributions, recursively aligning the each non-reference beam to aphysically closest beam towards the reference beam until aligning to thereference beam; and for the intensity distributions, recursivelyaligning the each non-reference beam to the physically closest beamtowards the reference beam comprises: aligning Beam_(i+1) to Beam_(i)and Beam_(i−1) to Beam_(i); aligning Beam_(i+2) to Beam_(i+1) andBeam_(i−2) to Beam_(i−1); repeating the alignments until each of thenon-reference beams is aligned to another beam towards the referencebeam; and propagating the alignments throughout the non-reference beamsto align each of the non-reference beams to the reference beam.
 2. Themethod of claim 1, wherein: the physical proximity among theintersections of the beams with the area corresponds to another physicalproximity among the light emitters on the LIDAR device.
 3. The method ofclaim 1, wherein: in association with the each beam, the each lightemitter and associated light detector are movable relative to the areaand configured to obtain an intensity reading at various time points toobtain the corresponding intensity distribution over the area.
 4. Themethod of claim 3, wherein: the each light emitter and associated lightdetector are configured to rotate relative to the vehicle as the vehiclemoves relative to the area; and in association with the each beam, theeach intensity distribution over the area comprises a plurality of theintensity readings caused by the rotation and the vehicle movement. 5.The method of claim 1, wherein: each of the intensity distributions isrepresentable in a histogram associating an light intensity reading anda frequency of occurrence of the light intensity reading in the area;the area comprises a marker corresponding to one or more intensity peaksin one or more histograms of the intensity distributions; and thereference intensity distribution is representative of the marker.
 6. Themethod of claim 1, further comprising: obtaining an intensity mappingthat maps any raw non-reference beam intensity reading to a calibratedintensity reading.
 7. A non-transitory computer-readable storage mediumcoupled to a processor and comprising instructions that, when executedby the processor, cause the processor to perform a method for lightintensity calibration, the method comprising: obtaining a plurality ofintensity distributions of reflected light from an area, wherein each ofthe intensity distributions is associated with a beam; determining areference intensity distribution from the plurality of intensitydistributions, wherein the reference intensity distribution isassociated with a reference beam, the plurality of intensitydistributions excluding the reference intensity distribution arenon-reference intensity distributions, and the non-reference intensitydistributions are each associated with a non-reference beam; andaligning each of the non-reference intensity distributions to thereference distribution by applying a Dynamic Time Warping (DTW)algorithm to calibrate the non-reference intensity distributions againstthe reference distribution, wherein: a LIDAR (Light Detection AndRanging) device comprises N rotatable light emitters and N detectorsmountable on a movable vehicle, the each light emitter emitting the eachbeam and the each light detector detecting the each corresponding lightintensity while in rotation as the vehicle moves relative to the area;the beams are Beams 1, 2, . . . i−1, i, i+1, i+2, . . . , and i+(N−i) inan order based on a physical proximity among intersections of the beamswith the area at the same time, i being any number from 1 to N; thereference beam is Beam_(i); aligning each of the non-reference intensitydistributions to the reference distribution comprises, for the intensitydistributions, recursively aligning the each non-reference beam to aphysically closest beam towards the reference beam until aligning to thereference beam; and for the intensity distributions, recursivelyaligning the each non-reference beam to the physically closest beamtowards the reference beam comprises: aligning Beam_(i+1) to Beam_(i)and Beam_(i−1) to Beam_(i); aligning Beam_(i+2) to Beam_(i+1) andBeam_(i−2) to Beam_(i−1); repeating the alignments until each of thenon-reference beams is aligned to another beam towards the referencebeam; and propagating the alignments throughout the non-reference beamsto align each of the non-reference beams to the reference beam.
 8. Thenon-transitory computer-readable storage medium of claim 7, wherein: thephysical proximity among the intersections of the beams with the areacorresponds to another physical proximity among the light emitters onthe LIDAR device.
 9. The non-transitory computer-readable storage mediumof claim 7, wherein: in association with the each beam, the each lightemitter and associated light detector are movable relative to the areaand configured to obtain an intensity reading at various time points toobtain the corresponding intensity distribution over the area.
 10. Thenon-transitory computer-readable storage medium of claim 9, wherein: theeach light emitter and associated light detector are configured torotate relative to the vehicle as the vehicle moves relative to thearea; and in association with the each beam, the each intensitydistribution over the area comprises a plurality of the intensityreadings caused by the rotation and the vehicle movement.
 11. Thenon-transitory computer-readable storage medium of claim 7, wherein:each of the intensity distributions is representable in a histogramassociating an light intensity reading and a frequency of occurrence ofthe light intensity reading in the area.
 12. The non-transitorycomputer-readable storage medium of claim 11, wherein: the areacomprises a marker corresponding to one or more intensity peaks in oneor more histograms of the intensity distributions; and the referenceintensity distribution is representative of the marker.
 13. Thenon-transitory computer-readable storage medium of claim 7, wherein themethod further comprises: obtaining an intensity mapping that maps anyraw non-reference beam intensity reading to a calibrated intensityreading.
 14. A system for light intensity calibration, comprising: Nrotatable light emitters and N detectors mountable on a movable vehicle,the each light emitter emitting a beam and the each light detectordetecting an reflection of the corresponding beam while in rotation asthe vehicle moves relative to an area; and a processor configured toreceive signals from the detectors and to: obtain a plurality ofintensity distributions of reflected light from the area, wherein eachof the intensity distributions is associated with the beam; determine areference intensity distribution from the plurality of intensitydistributions, wherein the reference intensity distribution isassociated with a reference beam, the plurality of intensitydistributions excluding the reference intensity distribution arenon-reference intensity distributions, and the non-reference intensitydistributions are each associated with a non-reference beam; and aligneach of the non-reference intensity distributions to the referencedistribution by applying a Dynamic Time Warping (DTW) algorithm tocalibrate the non-reference intensity distributions against thereference distribution, wherein: the beams are Beams 1, 2, . . . i−1, i,i+1, i+2, . . . , and i+(N−i) in an order based on a physical proximityamong intersections of the beams with the area at the same time, i beingany number from 1 to N; the reference beam is Beam_(i); aligning each ofthe non-reference intensity distributions to the reference distributioncomprises, for the intensity distributions, recursively aligning theeach non-reference beam to a physically closest beam towards thereference beam until aligning to the reference beam; and for theintensity distributions, recursively aligning the each non-referencebeam to the physically closest beam towards the reference beamcomprises: aligning Beam_(i+1) to Beam_(i) and Beam_(i−1) to Beam_(i);aligning Beam_(i+2) to Beam_(i+1) and Beam_(i−2) to Beam_(i−1);repeating the alignments until each of the non-reference beams isaligned to another beam towards the reference beam; and propagating thealignments throughout the non-reference beams to align each of thenon-reference beams to the reference beam.
 15. The system of claim 14,wherein: the physical proximity among the intersections of the beamswith the area corresponds to another physical proximity among the lightemitters.
 16. The system of claim 14, wherein: in association with theeach beam, the each light emitter and associated light detector aremovable relative to the area and configured to obtain an intensityreading at various time points to obtain the corresponding intensitydistribution over the area.
 17. The method of claim 16, wherein: theeach light emitter and associated light detector are configured torotate relative to the vehicle as the vehicle moves relative to thearea; and in association with the each beam, the each intensitydistribution over the area comprises a plurality of the intensityreadings caused by the rotation and the vehicle movement.
 18. The methodof claim 14, wherein: each of the intensity distributions isrepresentable in a histogram associating an light intensity reading anda frequency of occurrence of the light intensity reading in the area.19. The method of claim 18, wherein: the area comprises a markercorresponding to one or more intensity peaks in one or more histogramsof the intensity distributions; and the reference intensity distributionis representative of the marker.
 20. The method of claim 14, wherein theprocessor is further configured to: obtain an intensity mapping thatmaps any raw non-reference beam intensity reading to a calibratedintensity reading.